BABAR-PUB-10/031 
SLAC-PUB-14400 



Measurement of the mass and width of the I?si(2536)"'" meson 

J. P. Lees/ V. Poireau/ E. Prencipe,^ V. Tisserand/ J. Garra Tico,^ E. Grauges,^ M. Martinelli'^'',^ 
D. A. Milanes''^3 a. Palano"^3 M. Pappagallo''^3 G. Eigen," B. Stugu,"* L. Sun,"* D. N. Browii,^ L. T. Keitu, 

T. Schroeder,^ D. J. AsgeirssonJ C. HeartyJ 
I. Buzykaev,^ V. P. Druzhinin,^ V. B. Golubev,^ 
<^Vr..rr.c.r. 9 E. P. Solodov,^ K. Yu. Todyshev,^ 
° M. Mandelkern/o D. P. Stoker," 



Kerth,^ 

Yu. G. Kolomensky,^ G. Lynch,^ I. L. Osipcnkov,^ H. Koch,^ T. Schroeder,^ D. J. AsgeirssonJ C. H 
T. , 

E. A. Kravchenko,^ A. P. Onuchin,^ S. I. Serednyakov,^ Yu. I. Skovpen,^ E. P. Solodov,^ K. Yu. Todyshe-' 

. 10 



H. Atmacan," J. W. Gary," F. Liu," O. Long," G. M. Vitug," C. Campagnari, 

^ S 



S. 
M. Lo 

C. L. Lee, 26 M. ivium, rt. 

H. M. Lacker,29 T. Lueck,29 , , , , , 

H. B. Crawley,32 W. T. Meyer,32 S. Prell,32 E. L Rosenberg,32 A. E. Rubin,32 Gritsan,33 Z. J. Guo,33 

N. Arnaud,34 M. Davier,^^ D. Derkach,^^ J. Firmino da Costa,^^ G. Grosdidier,^^ F. Le Diberder,34 A. M. Lutz,^^ 
B. Malaescu,34 A. Perez,^^ p. Roudeau,^* M. H. Schune,^^ A. Stocchi,^^ L. Wang,34 q Wormser,34 D. J. Langc,^^ 

D. M. Wright,35 L Binghani,^^ C. A. Chavez,^^ J. p. Colenian,^^ J. R. Fry,^^ E. Gabathuler,^^ D. E. Hutchcroft,^^ 
D. J. Payne, C. Touramanis,^'' A. J. Bevan,^'^ F. Di Lodovico,^^ R. Sacco,^^ M. Sigamani,^'^ G. Cowan, 

S. Paramesvaran,38 A. C. Wren,^^ D. N. Brown,39 C. L. Davis,^^ A. G. Denig,40 M. Fritsch,^^ Gradl,^" 
A. Hafner,40 K. E. Alwyn,*! D. Bailey,^! R. J. Barlow,"! G. Jackson,"! G. D. Lafferty,"! R. Cenci,"2 B. Hamilton,"2 
A. Jawahery,"*2 D. A. Roberts,'*2 G. Simi,42 C. Dallapiccola,"^ E. Salvati,"^ R. Cowan,"" D. Dujmic,"" G. Sciolla,"" 
D. Lindcmann,"^ P. M. Patel,"^ S. H. Robertson,"^ M. Schram,"^ P. Biassom"^"6 A. Lazzaro°^"*^ V. Lombardo'',"^ 

F. Palombo''^"'^ S. Stracka°^"6 L. Cremaldi,"^ R. Godang,"^'0 R. Kroeger,"^ P. Sonnek,"^ D. J. Summers,"^ 
X. Nguyen,"^ P. Taras,"^ G. De Nardo°^"9 D. Monorchio''^"9 G. Onorato''^"9 C. Sciacca"\"9 G. Raven,^^ 

H. L. Snoek,50 C. P. Jessop,^! K. J. Knocpfcl,^! J. M. LoSecco,^! W. F. Wang,^! L. A. Corwin,52 K. Honscheid,52 
R. Kass,'52 Blount,53 Brau,^^ ^ Frey,^^ a. Kolb,^^ j^^ Rahmat,^^ q Sinev,^^ j^,^ Strom,^^ 

J. Strubc,53 E. Torrence,53 q Castclli"^^" E. Feltresi"^5" N. Gagliardi°^5" M. Margoni°\54 M. Morandin",'^" 
M. Posocco",^" M. Rotondo",^" F. Simonetto"^54 ^ Stroili°^^" E. Ben-Haim,'^^ jyj^ Bomben,^^ G. R. Bonneaud,^^ 

H. Briand,^^ G. Caldcrini,^^ J. Chauveau,^^ O. Hamon,^^ Ph. Leruste,^'^ G. Marchiori,^'^ J. Ocariz,^^ S. Sitt,'^^ 
M. Biasini'^^56 E. Manoni'^^56 A. Rossi'^^56 C. Angelini°^5'^ G. Batignani'^^57 Bettarini"^5'^ M. Carpinelli°^5^'[l 

G. Casarosa''^" A. Cervelli°^" F. Forti''^" M. A. GiorgP^57 a. Lusiani^^^T NerP^5^ E. Paoloni°^5^ 

G. Rizzo''^" J. J. Walsh*^," D. Lopes Pegna,^** C. Lu,^^ J. Olsen,^^ A. J. S. Smith,^^ A. V. Telnov,^* F. AnullP,^^ 
G. Cavoto",^^ R. Faccini"^59 F. Ferrarotto",^^ F. Ferroni''^59 M. Gaspero"^^9 L. Li Gioi°,59 M. A. Mazzoni'^,^^ 
G. Piredda",^^ C. Buenger,60 T. Hartmann,60 T. Leddig,^" H. Schr6der,60 R. Waldi,^" T. Adye,^! E. O. Olaiya,^! 
F. F. Wilson,^! S. Emery,62 q Hamel de Monchenault,62 q Vasseur,62 Yeche,*^2 x. Allen,^^ D. Aston,63 
D. J. Bard,63 Bartoldus,^^ j p Benitez,^^ (j Cartaro,^^ M. R. Convery,^^ j Dorfan,^^ q p Dubois-Felsmann,'^^ 

W. Dunwoodie,63 ^ q Field,63 M. Franco Sevilla,63 B. G. Fulsom,^^ a. M. Gabareen,^^ M. T. Graham,^^ 
P. Grenier,63 C. Hast,63 ^ ^ Innes,^^ M. H. Kelsey,63 H. Kim,63 p. Kim,63 M. L. Kocian,63 D. W. G. S. Leith,63 
P. Lewis,63 g Li 63 B Lindquist,63 g Luitz,63 V. Luth,63 H. L. Lynch,63 p, g MacFarlanc,63 p, Muller,63 
Neal,S3 g^ Nelson,63 C. P. 0'Grady,63 ofte,'^3 Pcrl,63 t. Pulliam,63 B. N. Ratcliff,'^3 g^ jj Robertson 
A. Roodman,63 a. A. Salnikov,63 y. Santoro,63 ^ jj. Schindler,63 J. Schwiening,'^3 a. Snyder,'^3 p,^ g^ 63 
M. K. Sullivan,63 S. Sun,63 K. Suzuki,63 j j^p Thompson,63 J. Va'vra,63 A. P. Wagner,63 M. Weaver,63 
W. J. Wisniewski,63 M. Wittgen,*^3 p, jj Wright,63 H. W. Wulsin,'^3 a. K. Yarritu,'53 C. C. Young,63 
V. Ziegler,63 X. R. Chen,^" W. Park,^" M. V. Purohit,^" R. M. Whitc^" J. R. Wilson,*^" A. Randle-Conde,^^ 
S. J. Sckula,65 M. Bcllis,^'^ P. R. Burchat,^^ T. S. Miyashita,^^ M. S. Alam,^^ J. A. Ernst,^^ N. Guttman,^^ 
A. Soffer,68 P. Lund,*^^ S. M. Spanicr,'^^ R. Eckmann,^" J. L. Ritchie,™ A. M. Ruland,™ C. J. SchiUing,™ 



2 



R. F. Schwitters,™ B. C. Wray7° J. M. Izen/i X. C. Lou/^ F. Bianchi°^^2 j)_ Ga^lba''^^2 PelliccionP\72 

L. LaIlceri''^'^3 Vitale°^'^3 jv^ Lopez-March/'' F. Martinez- Vidal/"* A. Oyanguren/'' H. Ahmcd/^ 
J. Albert/^ Sw. Banerjee/^ H. H. F. Choi/'^ k. Hamano/^ Q. J. King/^ R. Kowalcwski/^ M. J. Lcwczuk/^ 
C. Lindsay /5 I. M. Nugent/^ J. M. Roney/^ R. j. Sobie/'^ T. J. Gershon/^ p. F. Harrison,™ T. E. Latham,™ 
E. M. T. Puccio,™ H. R. Band,^^ S. Dasu,^^ Y. Pan,^^ R. Prepost,^^ C. O. Vuosalo,^^ and S. L. Wu^^ 

(The BaBAR Collaboration) 

' Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), 
Universite de Savoie, CNRS/IN2P3, F-74941 Annecy-Le-Vieux, France 
^ Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 
^INFN Sezione di Bari^ ; Dipartimento di Fisica, Universitd di Bari' , 1-70126 Bari, Italy 
'* University of Bergen, Institute of Physics, N-5007 Bergen, Norway 
^Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 
^ Ruhr Universitat Bochum, Institut fiir Experimentalphysik 1, D-44T80 Bochum, Germany 



7 



University of British Columbia, Vancouver, British Columbia, Canada V6T IZl 



8 I 



Brunei University, Uxbridge, Middlesex UBS 3PH, United Kingdom 
^Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 
"'University of California at Irvine, Irvine, California 92697, USA 
University of California at Riverside, Riverside, California 92521, USA 
^'''University of California at Santa Barbara, Santa Barbara, California 93106, USA 
University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 
'^California Institute of Technology, Pasadena, California 91125, USA 
'^University of Cincinnati, Cincinnati, Ohio 45221, USA 
'^University of Colorado, Boulder, Colorado 80309, USA 
'''Colorado State University, Fort Collins, Colorado 80523, USA 
'^ Technische Universitat Dortmund, Fakultdt Physik, D-44221 Dortmund, Germany 
''' Technische Universitat Dresden, Institut fiir Kern- und Teilchenphysik, D-0 1062 Dresden, Germany 
^" Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, F-91128 Palaiseau, France 
^' University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 
^^INFN Sezione di Ferrara"' ; Dipartimento di Fisica, Universitd di FerrarJ' , 1-44100 Ferrara, Italy 

^^INFN Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 
^'^INFN Sezione di Genova"^ ; Dipartimento di Fisica, Universitd di Genova^ , I-I6I46 Genova, Italy 
^''Indian Institute of Technology Guwahati, Guwahati, Assam, 781 039, India 
^^Harvard University, Cambridge, Massachusetts 02138, USA 
^'Harvey Mudd College, Glaremont, California 91711 
'''^ Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 
Humboldt- Universitat zu Berlin, Institut fiir Physik, Newtonstr. 15, D- 12489 Berlin, Germany 
^"imperial College London, London, SW7 2AZ, United Kingdom 
^' University of Iowa, Iowa City, Iowa 52242, USA 
^^lowa State University, Ames, Iowa 50011-3160, USA 
Johns Hopkins University, Baltimore, Maryland 21218, USA 
'^'^ Laboratoire de I'Accelerateur Lineaire, IN2P3/GNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 Orsay Cedex, France 
Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
'''University of Liverpool, Liverpool L69 7ZE, United Kingdom 
Queen Mary, University of London, London, El 4NS, United Kingdom 
University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 
'"''^University of Louisville, Louisville, Kentucky 40292, USA 
^"Johannes Gutenberg- Universitat Mainz, Institut fiir Kernphysik, D-55099 Mainz, Germany 
^' University of Manchester, Manchester M13 9PL, United Kingdom 
University of Maryland, College Park, Maryland 20742, USA 
University of Massachusetts, Amherst, Massachusetts 01003, USA 
''^'^ Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 
^"McGill University, Montreal, Quebec, Canada H3A 2T8 
INFN Sezione di Milano""; Dipartimento di Fisica, Universitd di Milano'', 1-20133 Milano, Italy 
'^''University of Mississippi, University, Mississippi 38677, USA 
Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 3J7 
^''INFN Sezione di Napoli"' ; Dipartimento di Scienze Fisiche, 
Universitd di Napoli Federico if', 1-80126 Napoli, Italy 
^" NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

"' University of Notre Dame, Notre Dame, Indiana 46556, USA 
^^Ohio State University, Columbus, Ohio 43210, USA 
University of Oregon, Eugene, Oregon 97403, USA 



3 



^^INFN Sezione di Padova"' ; Dipartimento di Fisica, Universita di PadovJ' , 1-35131 Padova, Italy 
''^ Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Marie C'urie-Paris6, 
Universite Denis Diderot- Paris!, F-75252 Paris, France 
^^INFN Sezione di Perugia"" ; Dipartimento di Fisica, Universita di Perugia^ , 1-06100 Perugia, Italy 
INFN Sezione di Pisa"" ; Dipartimento di Fisica, 
Universita di Piso!' ; Scuola Normale Superiore di Pisa'^ , 1-56127 Pisa, Italy 
'^^ Princeton University, Princeton, New Jersey 08544, USA 
^^INFN Sezione di Roma"" ; Dipartimento di Fisica, 
Universita di Roma La Sapienza^ , 1-00185 Roma, Italy 
Universitdt Rostock, D-18051 Rostock, Germany 
''^Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 
'^^CEA, Irfu, SPP, Centre de Saclay, F-91191 Gif-sur-Yvette, France 
^^SLAC National Accelerator Laboratory, Stanford, California 94309 USA 
^■^ University of South Carolina, Columbia, South Carolina 29208, USA 
^^Southern Methodist University, Dallas, Texas 75275, USA 
'"'Stanford University, Stanford, California 94305-4060, USA 
''''State University of New York, Albany, New York 12222, USA 
''^Tel Aviv University, School of Physics and Astronomy, Tel Aviv, 69978, Israel 
""University of Tennessee, Knoxville, Tennessee 37996, USA 
'"'University of Texas at Austin, Austin, Texas 78712, USA 
'" University of Texas at Dallas, Richardson, Texas 75083, USA 
'INFN Sezione di Torino"" ; Dipartimento di Fisica Sperimentale, Universita di Torino'' , 1-10125 Torino, Italy 
'''' INFN Sezione di Trieste" ; Dipartimento di Fisica, Universita di Trieste'' , 1-34127 Trieste, Italy 
^'^IFIC, Universitat de Valencia- CSIC, E-46071 Valencia, Spam 
'"'University of Victoria, Victoria, British Columbia, Canada V8W 3P6 
Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 
''''University of Wisconsin, Madison, Wisconsin 53706, USA 
(Dated: January 14, 2013) 

The decay width and mass of the _Dsi(2536)^ meson are measured via the decay channel 
D^^ — > D*'^ K1 using 385 fb^^ of data recorded with the BABAR detector in the vicinity of 
the T(46') resonance at the PEP-II asymmetric-energy electron-positron collider. The result for 
the decay width is r(Z)+) = 0.92 ± 0.03 (stat.) ± 0.04 (syst.) MeV. For the mass, a value of 
m(L»+J = 2535.08 ± 0.01 (stat.) ± 0.15 (syst.) MeV/ is obtained. The mass difference between the 
L>+i and the D'+ is measured to be m(D+i) - m(D*+) = 524.83 ± 0.01 (stat.) ± 0.04 (syst.) MeV/c^, 
representing a significant improvement compared to the current world average. The unnatural 
spin-parity assignment for the Df^ meson is confirmed. 



PACS numbers: 14.40.Lb, 13.25.Ft, 13.66.Bc 
I. INTRODUCTION 

The theoretical description of mesons is problem- 
atic because, unlike D mesons, the masses and widths 
of the D;q{2317)+ and Dsi{2460)+ states [H-Q are not 
in agreement with potential model calculations based on 
HQET Theoretical explanations for the discrepancy 
invoke D^*^K molecul es B, c hiral partners 0, [id , uni- 
tarized chiral models ll ll , tetraquarks [Ij, , and 
lattice calculations Il5lll6i. but a satisfactory description 
is still lacking (see [13, [l3 fo'" more details). Improved 
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measurements of the -Dj^ meson parameters can lead to 
a better understanding of these states. 

In this analysis a precise measurement of the 
Z?si(2536)"'" mass and decay width is performed based 
on a high statistics data sample [l^. The Dsi(2536)"'' 
meson, referred to as the D^^ below, was first seen in 
cc-continuum reactions [20j . and more recently in B 
decays. The current world average mass value pub- 
lished by the Particle Data Group is based on measure- 
ments with large statistical and systematic uncertainties: 
2535.29 ± 0.20 MeV/c^ [2l|; the mass difference between 
the Df^ and the D*~^ meson has been measured to be 
525.04 ± 0.22 MeV/c^ [U. An upper limit on the decay 
width (r < 2.3 MeV at 90% confidence level), and a mea- 
surement of the spin-parity of the D'^i meson ( = 1+), 
have been obtained, but based on low-statistics data sam- 
ples only [2l| - [23j . The mixing between the D'^^ meson 
and the other = l"*" state Dsi(2460)"'' was investi- 
gated in Ref. (H]. 

This analysis is based on a data sample corresponding 
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to an integrated luminosity of 349 fb~^ recorded at the 
T{4S) resonance and 36 fb~"^ recorded 40 MeV below that 
resonance with the BaBAR detector at the asymmetric- 
energy e+e" collider PEP-II at the SLAC National Ac- 
celerator Laboratory. In this analysis, D'^^ mesons are 
reconstructed from cc continuum events in the D*~^K'^ 
channel; those originating from B decays arc rejected. 

The BaBAR detector is described briefly in Sec. |lll 
The principal criteria used in the reconstruction of 
the D*^Kg mass spectrum and the selection of D^^- 
candidates are discussed in Sec. IIIII The relevant Monte 
Carlo (MC) simulations are described in Sec. llVi while 
the detector resolution parametrization is considered in 
Sec. |Vl Measurements of the mass and total width for 
the D^^ state are obtained from a fit to the D*^K'^ in- 
variant mass distribution as discussed in Sec. IVIl Decay 
angle distributions are studied in Scc. lVIll where the im- 
plications for the spin-parity of the Df^ state are also dis- 
cussed. Sources of systematic uncertainty are described 
in Sec. IVIIIl and the results of the analysis are summa- 
rized in Sees. llXl andlXl 



II. THE BABAR DETECTOR 

BaBAR detector is described in detail else- 
25| . Charged particles are detected, and their 



The 
where 

momenta measured, with a combination of five layers of 
double-sided silicon microstrip detectors (SVT) and a 40- 
layer cylindrical drift chamber (DCH), both coaxial with 
the cryostat of a superconducting solenoidal magnet that 
produces a magnetic field of 1.5 T. Charged particle iden- 
tification is achieved by measurements of the energy loss 
dE I dx in the tracking devices and with an internally re- 
fiecting, ring-imaging Cherenkov detector. The energy of 
photons and electrons is measured with a CsI(Tl) electro- 
magnetic calorimeter, covering 90% of the Att solid angle 
in the T(4S') rest frame. The instrumented flux return 
of the magnetic field is used to identify muons and A'°'s. 



III. SELECTION AND RECONSTRUCTION OF 
EVENTS 

The Df^ is reconstructed via its decay mode D*^K^, 



with K2 



■K^TT and D* 



D"tt+. The D° is 



reconstructed through two decay modes, K^t:^ and 
A'^TT+TT+TT^ , which will be labeled KAt: and A'Gtt, re- 
spectively, in the following. To improve the mass resolu- 
tion, the mass difference Am{D'^^) = m{Df^)—m{D*'^) — 
m{K'^) is examined. 

Events are selected by requiring at least five charged 
tracks, at least one of which is identified as a charged 
kaon. Also, at least one neutral kaon candidate is re- 
quired. Each track must approach the nominal e~^e~ 
interaction point to within 1.5 cm in the transverse di- 
rection, and to within 10 cm in the longitudinal (beam) 
direction. Kaon candidates are identified using the nor- 



malized kaon, pion and proton likelihood values {Lk, 
and Lp) obtained from the particle identification system, 
by requiring Lk /{Lk + L^^) > 0.50 and Lk/{Lk + Lp) > 
0.018. Furthermore, the track must be inconsistent with 
the electron hypothesis or have a momentum less than 
0.4GeV/c. Tracks that fulfiU Lk/{Lk + L^) < 0.98 and 
Lp/{Lp + Ltt) < 0.98 are selected as pions. 

Candidates for the decay arc formed by selecting 
all K^TT^ pairs (K~tt'^tt^tt~ combinations in the second 
mode) that have an invariant mass within ±100 MeV/c^ 
of the nominal mass [2lj. Candidates for the D*~^ de- 
cay are formed by adding a 7r+ to the such that 
the mass difference between and D° is less than 

170 MeV/ c^. A AT^ candidate consists of a 7r+ tt^ pair 
with invariant mass within ±25 MeV/c^ of the nominal 
mass [2l| . A kinematic fit is applied to the complete de- 
cay chain, constraining the D'^^ candidate vertex to be 
consistent with the e'^e" interaction region. Mass con- 
straints are not applied to intermediate states. Those D'^^ 
candidates with a fit probability greater than 0.1% 
are retained. To suppress combinatorial background and 
events from _B-decays, we require the momentum p* {D'^^ ) 
of the D^j^ in the T(4S') centcr-of-mass (CM) frame to ex- 
ceed 2.7GeV/c. 

The A'tt and Ktttttt mass spectra for accepted D° can- 
didates, shown in Figs, [ija) and[IJd), are fitted with a 
signal function consisting of a sum of two Gaussians with 
a common mean value, and a linear background function. 
The width of the signal regions for Z3*+ and AT" can- 
didates is defined as twice the full width at half maximum 
(FWHM) of the signal line shapes. A signal window of 
±18 (14) MeV/c^ for the A:47r (A:67r) mode around the 
mean mass of 1863.5 (1863.5) MeV/ obtained from the 
fit is used to select D'^ candidates. For these candidates, 
the ZJ^Tr"*" — D° mass difference distributions are shown in 
Figs-Htb) and[Tfe). These are fitted with the sum of a rel- 
ativistic Brcit-Wigner signal function and a background 
function consisting of a polynomial times an exponential 
function. A Z3*+ signal region of ±1 McV/c^ around the 
fitted mean value of 145.4 MeV/ is chosen for both de- 
cay modes. To further reduce the background, the angle 
between the flight direction of the K° candidate and the 
line connecting the e~^e~ interaction point and the A'^ 
decay vertex is required to be less than 0.15 radians. For 
candidates passing these selection criteria, the ATg candi- 
date invariant mass distributions (Figs.^c) andfljf)) are 
fitted with the sum of a signal function, consisting of the 
sum of two Gaussians, and a linear background function. 
A signal window of ±6MeV/c^ around the fitted mean 
mass of 497.2 MeV/c^ is selected for both decay modes. 

In the case of an event with multiple candidates, the 
candidate with the best fit probability is chosen. The 
Am{D'^^) candidate spectra after all selection criteria are 
shown in Figs. [^Ka) and [2Kb). The fits to these spectra 
use a Double-Gaussian signal function and a linear back- 
ground function. Note that for this preliminary fit the 
intrinsic width and the resolution are not taken into ac- 
count. The FWHM values obtained are (2.2 ± 0.1) MeV 
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and (2.0 ± 0.1) McV, respectively, with corresponding 
signal yields of about 3500 and 4000 entries. 



IV. MONTE CARLO SIMULATION AND 
COMPARISON WITH DATA 

Monte Carlo events are generated for Df^ — !> 
D*+Kl,D*+ D°Tr+, with D" K-tt+ and 
K^TT^TT^n^ , by EvtGen [2^. The detector response is 
simulated using the GEANT4 package. For each D" 
decay mode, and for each of the corresponding D~i de- 
cays, 776000 events are generated. The D^^ line shape is 
generated using a non-relativistic Breit-Wigner function 
having central value m{D^^)gen = 2535.35 MeV/c^ and 
intrinsic width T{D'^^)gf,n = 1 MeV (this sample is labeled 
Fi in the following). The range of generated D^^ masses 
is restricted to values between m{Df^)gen — lOMeV/c^ 
and m(D^j^)ge„ + 15 MeV/c^. The masses of the daughter 
particles are taken from Ref. [2l| . 

In order to test the mass resolution model, a second 
set of MC samples with 381000 events for each de- 
cay mode is generated using a Breit-Wigner width of 
^{Dfi)gen = 2 MeV (F2 sample). In addition to these 
signal MC samples, separate D° and K'^ samples are cre- 
ated from data and generic cc MC simulations without 
requiring a D*'^ or D^^. They are used mainly for reso- 
lution studies. 

The MC and data arc in good agreement for the trans- 
verse momentum distributions of pions, kaons, D and 
D* mesons, and for the number of SVT coordinates of 
pions and kaons. The agreement is worse for the number 
of DCH coordinates, where the data show systematically 
fewer coordinates than the MC, giving rise to a resolu- 
tion that is about 10% smaller in the MC than in data. 
This is illustrated in Fig. 121 which shows the p*{K'^) and 
p*{D^) dependence of the ratio between the FWHM of 
the resolution functions in cc MC and data, where p* 
is the momentum in the CM frame. This effect will be 
discussed further in Sec. IVIIII There is also disagree- 
ment between the number of D^^ signal entries in MC 
and data as a function of p*[D^^) (Fig. |4]). This effect 
will be addressed in Sees. IVl and I Villi 



V. RESOLUTION MODEL 

The resolution model is derived from the D'^-^ signal 
MC by studying the difference Am^es between the re- 
constructed and generated mass values. The Multi- 
Gaussian ansatz 

G'(Am^es) = / -^e ^ dcr (1) 

is found to accurately model the mass resolution spec- 
tra. This represents a superposition of Gaussian distri- 
butions with the same mean value Am^eso but variable 
width (T, starting from minimum width ctq a-nd increasing 



TABLE I: Reconstructed values for A?n(Z)+i)o and T{D+^) 
(fit to MC sample Fi). The resolution model used is derived 
from MC sample Fi. 





Am(D+)o/MeV/c2 


r(D+)/MeV 




27.737 ± 0.003 


1.001 ±0.005 


{KQtv) 


27.734 ± 0.003 


0.991 ± 0.006 



to maximum width ruo • The FWHM of the distribution 
is numerically calculable once (Tq and r are known. The 
mass resolution for the different particles depends on the 
CM momentum p*{D'^^). Therefore, the parameter (Tq of 
Eq. (|T]) is obtained as a function of p*{D^^). 

Figuresini^a) and|5Kb) show Am^es distributions for the 
full p*{D'l^) range. From these plots the value of the pa- 
rameter r is determined to be 4.78 ± 0.04 and 5.20 ± 0.05 
for the KAtt and KQii modes, respectively. Events 
are divided into 30 p*{D'^^) intervals from 2.7GeV/c to 
4.7GeV/c and the fit repeated for each interval, result- 
ing in p* (Z? )-dependent ctq values (Figs.lHJa) andjBJb)). 
The corresponding p*(D^j^)-dependent FWHM of the res- 
olution functions is shown in Figs.jTla) andjZl^b). 

In order to validate this resolution model, the p*{D^T^)- 
dependent resolution function with the corresponding pa- 
rameters ctq and r is convolved with a non-relativistic 
Breit-Wigner function and fitted to the Am(Df^) sig- 
nal MC distribution (MC sample Fi). The results are 
shown in Figs.jSJa) andjSJb), and the reconstructed val- 
ues for mean A'm{D'^^)Q and width T{D^^) are listed in 
Table HI The corresponding generated values for both 
decay modes are Am{Df^^)gen = 27.744 MeV/c^ for the 
mean and r{Dj^)gen = 1.000 MeV for the width. The 
small deviations between generated and reconstructed 
values are discussed in Sec. IVIlH 



VI. FIT TO THE D'K^ MASS SPECTRUM 

For the final fit to the D*K'^ mass spectra, as repre- 
sented by the Am{D^^) distributions of Figs. [Hand |9l the 
signal function consists of a relativistic Breit-Wigner line 
shape numerically convolved with the p* (I?^x )-dcpendent 
resolution function (Eq. ([T])). A linear function is used 
to describe the background. 

The relativistic Breit-Wigner function used takes the 
form 

VPi.mo / \m ) {ml- m?Y + ^l^ml ' 

where mo is an abbreviation for Am{D^^)[) and m stands 
for Am{D'^^). The variable pi,m is the momentum of the 
D*'^ in the rest frame of the D'^^ resonance candidate, 
which has mass to, and Pi,mo is the value for m = toq. 
The respective Blatt-Weisskopf barrier factors F]^{pi „i) 
for orbital angular momentum L between the D*~^ and 
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FIG. 1: (a,d) Z?" candidate invariant mass distributions; (b,e) Difference between the Z)*^ and D° candidate invariant masses. 
(c,f) Ks candidate invariant mass distributions; (a-c) Kin mode; (d-f) K6n mode. Signal regions are indicated by tlie vertical 
lines. The signal and background line shapes fitted to the mass distributions are described in tlie text. 



K2 are 



FoiPi,m) 

F2{pi,m) 



= 1, 



(3) 
(4) 

(5) 
(6) 

is defined as in Ref. psj . Th.e mass-dependent widtli is 
given by 



where 



v/l + (flpi.rnoF 

^9 + SiRpunY + {RpunY 



i? = 1.5 (GeV/c)-i 




FL{p2.mf 



(7) 



with T{D^i) the total intrinsic width of the D^i res- 
onance. This relation takes into account the D^^ 
D*+K° and the D+i D*°K+ decay modes, with the 
corresponding branching fractions Bi and B2, respec- 
tively: 



B, = 



p2L + l 



2L+1 , 2L+1 ■ 
Pl,mo ^ P2,ma 



(8) 



Since the D'^^ mass lies close to threshold for both de- 
cay modes, the mass values of the decay particles make 
a significant difference. The momenta P2,m and P2,ma 
correspond to pi,m and pi.moi respectively, but are cal- 
culated for the D*^K^ decay mode. 

It is assumed that the D'^-^ has spin-parity = 1+ and 
from parity conservation that the orbital angular momen- 
tum L is either or 2. The S'-wave usually dominates 
in 1+ decays, so L = is chosen for the main fit and 
an additional L ~ 2 contribution is used to estimate a 
systematic uncertainty. Further discussion on the J and 
L values is presented in Sec. IVIII 

The fit to the Am{D+j^) = m(D+ ) - 7n(D*) - m{K°) 
mass difference spectrum in data (Fig. ^ yields mean 
mass differences 

Am(Z3+i)o = 27.231 ± 0.020 MeV/c^ (KAn), 
Am(D+)o = 27.205 ± 0.018 MeV/c^ (KGir), 

and total width values 

rp+i) = 1.000 ± 0.049 MeV (KAir), 
rp+i) = 0.941 ±0.045 MeV {K6tt). 

The fitted values for the two decay modes agree 
within the statistical errors. The signal yield is 3704 ±71 
for KiTT and 4334 ± 78 for K6n. 
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FIG. 2: Am(D+ ) = m(D+ ) - m(D*+) - m(A'g) invariant 
mass distributions in data after applying all selection criteria 
for the (a) KAn and (b) KGir mode. A Double- Gaussian signal 
function and a linear background function are used to describe 
the data in a preliminary fit. 



VII. ANGULAR DISTRIBUTION 

The assigned spin-parity = 1+ of the D^-^ is 
based on studies with small data samples (less than 200 
reconstructed events) [13, [111. There, fits of an angular 
distribution corresponding to unnatural spin-parity 
(1+, 2~, . . .) yielded the highest confidence level. In this 
analysis clean signals with a total number of about 8000 
reconstructed Dj"]^ -candidates are available, making a 
detailed study possible. 



decay angle. Since in this analysis the origin 
of the D^^ is not known, the decay angle 6' between the 
_D" momentum vector in the D*'^ CM system and the 
momentum vector in the CM system (Fig. [TUk ) 
is used for the analysis. The resulting angular dis- 
tribution dN {D'^^) / d cos 6' is infiuenced by the spin of 
the D'^^. The expected distributions for different D'^^ 
spin-parity values are calculated using the helicity for- 
malism j29l - [3l| and are listed in Table HIl 

The data are corrected for the detection efficiency and 
divided into 20 bins of cos 9'. The signal entries for the 
cos 9' bins are obtained from separate fits to the data 
with the mass and decay width of the Df^ fixed to the 



FIG. 3: ^'-dependence of the ratio between the FWHM of the 
resolution functions from cc-MC and data, (a) tt+tt". 
(b) D'' — >■ if ^7r^7r"'"7r~ . The solid line shows the fitted mean 
ratio with a value of 0.9. 

values reported in Sec. |Vll The dN{D+^)/d cos 9' distri- 
bution shown in Fig. [TT] is the combined result from the 
KAtt and A'Gtt samples. 

Comparison with the theoretical distributions shows 
a clear preference for the unnatural spin-parity values 
= l+,2^,3+..., confirming the earlier results [l^, 
[23} . The signal function for these values is 

Ii9') = a{sm^9' + Pcos^9'), (9) 

where /3 — |ylooP/|^ioP and a is a constant. The 
helicity amplitudes \Aoq\ and |^io| correspond to the 
helicities and ±1, respectively. 
The lowest value = 1+ is the most probable one: 
assuming 1"*" implies / = 1 (orbital momentum between 
the light and heavy quark), while the higher J values 
demand I > 2; such mesons are expected to be highly 
suppressed in e~^e~ production [32|. The coefficient 
/3 is 1 in the case of a pure S-wavc decay of D'^-^ into 
Z?*^/^^. thus yielding a flat distribution in disagreement 
with data. The results reported in Table HIl for /3 clearly 
indicate a D-wave contribution. Based on the results 
for /?, the ratio of the helicity amplitudes is determined 
to be |^io|/|Aoo| = 2.09 ± 0.09 for the combined Kin 
and K6tt samples, and 2.09 ± 0.14 and 2.04 ± 0.13 
for the individual samples, respectively. The squared 
ratio of the amplitudes is |AioP/|AooP = 4.35 ± 0.38 
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FIG. 4; p*(Z)^j^)-dependence of the D^j^ signal yield for data 
(open squares) and reconstructed MC (solid points) for the 
(a) K4:TT and (b) KGn decay modes. 
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FIG. 5: Fit of the resolution function (Eq. {T])) to Arrires with 
the r and ao parameters free to vary for the (a) K4-k and (b) 
K6tt decay modes. 
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2 = 3.6 ±0.3 ±0.1 



with 

SI. 



the BeUe resuh 



Several effects that might affect the resuhs of the an- 
gular analysis have been studied. 



D+i decay angle. The dN{Dj^)/d cos 9 distribu- 
tion is also studied, where 9 is the decay angle be- 
tween the Z3*+ momentum vector in the Df^ CM system 
and the D'^^ momentum vector in the e^e~ CM system 
(Fig. llOb). The combined efficiency-corrected cos6' spec- 
trum is shown in Fig. [121 The results in this figure indi- 
cate that the Df^ decay to D*~^K^ is not purely 5-wave. 
Were this decay purely S'-wave, the distribution would 
be flat. The cos 61 distribution, assuming = 1+, is 

I{9) = a((l+poo)|^io|' + (l-poo)|^ooP 

+ (l-3poo)(|^io|'-|^oo|')cos2 0) (10) 

where poo gives the probability that the ZJ^j helicity is 
zero. 

Results from a fit of both a constant and a distribution 
proportional to 1 + 1 cos^ 9 (based on Eq. ([TU)) ) are given 
in Table Hill Using the value of t from the cos 6* fit, the 
result for |AioP/|AooP from the cos 9' fit, and the coeffi- 
cients from Eq. ^TO\i . we determine poo = 0.48 ± 0.03 for 
the combined KAtt and K6tt samples, and 0.44±0.04 and 
0.53 ±0.04 for the individual samples, consistent with the 
Belle resuh poo = 0.490 ± 0.012 0. 



Test for non-flat efficiency. The formalism used 
for the calculation of I{9') assumes a flat acceptance in 
cos 9. In this study, the efficiency decreases a few percent 
for values of cos 9 > 0. In order to assess the impact of 
this effect, all D'^^ candidates with cos 9 > are removed 
from the data sample. The results for /? from fits to the 
reduced cos 9' spectra are consistent with the nominal 
results, ruling out an observable effect due to non-flat 
efficiency. 

Test for possible interference. Possible interfer- 
ence with unreconstructed recoil particle(s) X in the de- 
cay chain e+e~ ^ti-^ considered. The effect of 
interference is expected to depend on the flight direction 
of the D^^ . Therefore the data are divided into four sub- 
samples based on their cos 6^ value, where 9d is the flight 
angle of the Df^ relative to the beam axis (calculated in 
the e+e" CM system). For each of these reduced data 
samples, the fit to the cos 9' distribution is repeated. The 
values obtained for the parameter /3 are fully consistent 
within errors with each other and with the nominal value 
(full data sample) , ruling out a significant interference ef- 
fect. The same consistency between results is found in 
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fits to cos 0. 



VIII. SYSTEMATIC STUDIES 

The investigated sources of systematic uncertainty can 
be separated into three main categories: uncertainties 
arising from the resolution model, fit procedure, and re- 
construction. The uncertainties are defined by taking the 
differences Aajti and Ar between the standard result for 
the mass difference Am(D+i)o and width T{D+j^) given 
in Sec. I VII and the result obtained with the correspon- 
dent modification. A summary of the results is listed 
in Tabic IIVI If not otherwise stated, the momentum- 
dependent resolution model and the rclativistic Breit- 
Wigncr signal function combined with a first order poly- 
nomial for background parametrization from the stan- 
dard fit are used. Deviations smaller than 0.5keV/c^ for 
Aato and smaller than 0.5 keV for Ar are considered as 
negligible. 



(Fig. [3]), yielding an overestimated decay width from the 
fits to data. The effect of this is quantitatively studied by 
increasing the width of the resolution function by 10%. 
The repeated fits yield no significant deviations for the 
mass difference, but a 51 (45) keV smaller decay width. 
The nominal decay width values obtained from the fits in 
Sec. IVII are thus corrected by these values, yielding val- 
ues of T{Df^) = 0.949 (0.896) McV for the KAn {K6tt) 
mode. 

To estimate the corresponding systematic uncertainty, 
the resolution function modification is varied within (10± 
4)% to take a possible p* dependence into account (this 
value is derived from Fig. ^a), which shows the largest 
variation in p*). There are no effects on Am{D'^^)o, and a 
deviation of keV for T{Df-^) is observed in both decay 
modes, compared to the corrected results from above. As 
a conservative estimate, the larger deviation is used as a 
two-sided uncertainty, providing the largest systematic 
error for the decay width. 



A. Resolution model uncertainties 

General comparison between MC and Data. 

The D'^ and test samples (see Sec. IIV|) demonstrate 
that the mass resolution is underestimated by 10% in MC 



Further validation of the resolution model. To 

further validate the procedure used to obtain the res- 
olution model, the results from fits to the Fi and F2 
MC samples are compared. The derived resolution func- 
tion parameters are in good agreement between the two 
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TABLE II: List of spin-parity values for the _Dj^ and the corresponding decay angle distributions for the D*^ . Under the 
assumption of a strong decay, 0^ is forbidden. The last three columns show the y^^ /NDF of the fits to the cos ^'-distribution 
for efficiency-corrected data, with NDF being the number of degrees of freedom. 

dN{D+^)/d cose' xVNDF{K4tv) /NDF{K6-!v) xV^^^J" (combined data) 

0^*^ forbidden ^ ^ ~ 

0" acos^e' 2142.7/19 2440.8/19 4578.0/19 

l",2+,3",... asin^e' 103.2/19 108.8/19 190.9/19 

1+, 2", 3+,... (5*- wave only) const 392.1/19 425.1/19 802.5/19 

l+,2-,3+,... L»-wave) a(sin2 6*' -(- ^ cos^ 5') 24.9/18 9.5/18 14.7/18 

= 0.23 ± 0.03) (/? = 0.24 ± 0.03) (/3 = 0.23 ± 0.02) 



TABLE III: /NDF values of the fits to the cos ^-distribution for efficiency-corrected data, with NDF being the number of 
degrees of freedom. 





dN{D+^)/d cose 


X^/NDF{KA-k) 


X^/NDF{K6tt) 


X^ /NDF (combined data) 


pure 5- wave 


constant 


19.0/19 


55.5/19 


57.0/19 


S- and D-wave 


a(l + tcos^ e) 


12.0/18 


27.3/18 


25.2/18 






{t = -0.15 ±0.05) 


{t = -0.27 ±0.05) 
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FIG. 8: Fit of a non-relativistic Breit-Wigner convolved with 
the resolution function to the _DJ^ candidate mass difi'erence 
spectra in the Fi MC sample for the (a) i('47r and (b) KQ-r 
decay modes. 



2.018 ± 0.022 McV for the K&i: mode, are in good agree- 
ment with the generated values. Similarly, when the res- 
olution function from the T2 sample is used to deter- 
mine the width for the Fi sample, values of T{D^^) = 
1.003 ± 0.005 MeV and 0.999 ± 0.006 MeV, respectively, 
are obtained, in agreement with the generated values. 

As a conservative estimate, the small deviations found 
during the validation procedure for the resolution model 
using the Fi sample in Sec. |V] arc used as systematic un- 
certainties: AAm = -7 (-10)keV/c2; Ar = +1 (-9)keV 
for KAt: {K6tt). 

Alternative resolution models. Using the resolu- 
tion model obtained from the F2 MC sample, a fit to 
data yields uncertainties Aatk < 0.5 (< 0.5)kcV/c^ and 
Ar = +1 (+12)keV for KAtt {K67t). 

Instead of the momentum-dependent resolution model 
of the standard analysis, an alternative model has been 
tested, based on the comparison of MC and data dis- 
tributions that show disagreement, such as the p*{D^^) 
dependence of the Df^ yield. By dividing the MC and 
data spectra from Fig. 21 a correction function is derived. 
MC data are modified with this function such that the 
two distributions in Fig. H] coincide. From these corrected 
MC, a new resolution model is derived. The results for 
Am(_D^j)o and F(Z?^j^) in data agree within the error with 
the momentum-dependent treatment (systematic uncer- 
tainties Aa™ < 0.5 (< 0.5)keV/c2, Ar = -2 (±l)kcV 
for KAn {K6tt)). 

The larger deviations listed above are reported as the 
systematic uncertainties associated with the use of alter- 
native resolution models. 



samples. The widths of the reconstructed Am{Df-^) Parameters of the p*(£)^)-dependent reso- 

distributions determined from fitting the F2 samples, lution model. The parameter r of the p*{D'^^)- 
r(-Dji) = 2.004 ± 0.016 MeV for the KAtt mode and dependent resolution model is modified within its er- 
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FIG. 9: Fit of a relativistic Breit-Wigner convolved with the 
resolution function to the Df^ candidate mass difference spec- 
tra in data, for the (a) K4-k and (b) K6n modes. The dotted 
line indicates the background line shape. The upper parts of 
the figures show the normalized fit residuals. 



FIG. 11: Efficiency-corrected signal yield as function of cos 9' 
in data. The following models are fitted to the distribution: 
a(sin^ 6' + 13 cos^ 0') (solid line); a constant (dash-dotted line); 
a cos^ 9' (dashed line) ; a siv? 9' (dotted line) . 



ror leading to negligible deviations in /^m{D^^)Q and 
±6 (±7)keV in T{D+^) for KAi: {K&tt). A different 
parametrization of the cro(p*(^^L))-dependence (second 
order polynomial) results in a negligible deviation for 
/\m{D+^)o and -3 (-2) keV for r(L>+ ). 

B. Fit Procedure Uncertainties 

Breit-Wigner line shape. In the standard fit, a 
pure 5-wave decay of the D^^ to D*^ is assumed, us- 
ing a Breit-Wigner line shape corresponding to L = 0. 
To estimate a systematic uncertainty, a combination of 
an /S-wave and a Z?-wave Breit-Wigner is used instead. 
Relative contributions of 72% and 28% are used, based on 
a decay angle analysis of the D^^ by the Belle Collabora- 




200 7 - 
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r ' 
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FIG. 12: Efficiency-corrected signal yield as function of cos^ 
in data. The following models are fitted to the distribution: 
constant (dotted line); a(l + t cos^ 9) (solid line). 
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tion [i^l • Using the modified signal lineshape, uncertain- 
ties of -9 (-8) keV/c^ in Am{D+^)o and -2 (-3) keV in 
r(D^j^) are derived, compared with, the standard results. 

As an additional check, the value of R (Eq. ([H])) is set 
to 2.0 (GeV/c)-i. No effect on Ato(D+ )o and r(D+ ) is 
observed. 

The effect of neglecting £)*°/v+ decays (Sec. |VT| is 
studied by setting Bi — 1 and B2 ~ 0. The resulting un- 
certainties arc negligible for both Am(D^j^)o and r(_D^]^). 

Numerical precision of convolution. The inte- 
gration range and number of steps in the numerical con- 
volution of the signal line shape and resolution function 
(Sec. IVip are varied, resulting in a negligible deviation 
both for the mass and the width. 

Mass window. The mass window for Am{D'^^) is 
enlarged, resulting in no significant change for Am(£)^j^)o 
and a difference for the width of Ar = +9 (+3) keV. 

Background parametrization. For background 
parametrization, a power law function proportional 
to m" is used instead of a linear function, leaving 
Am(Z)^j)o unaffected but yielding Ar = —5 (— 7)keV 
for KAtt {K6tt). 



C. Reconstruction Uncertainties 

Tracking region material. Uncertainties in the 
mass may arise from uncertainties in the energy- 
loss correction in charged particle tracking. Studies of A 
and K'^ decays suggest that the correction might be un- 
derestimated dl]. Two possibilities are considered, one 
with the SVT material density increased by 20% and the 
other with the tracking region material density (SVT, 
DCH) increased by 10%, as was investigated in detail in 
Rcfs. [l,!!^. The deviations indicate that the fit results 
for the mass might be underestimated. The larger val- 
ues from the two studies (Aato = +21 (+13) kcV/c^ and 
Ar = +14 ( — 15)kcV for KAtt (KGtt)) are chosen as a 
two-sided systematic uncertainty. 

SVT alignment. Slight possible deviations in the 
alignment of SVT components may affect the measure- 
ment of angles between tracks and thus the mass mea- 
surement. This is studied by applying small distortions 
to the SVT alignment in simulated data, comprising 
general changes between different run periods and ra- 
dial shifts. Results are Aajk = ±6 (±7)keV/c^ and 
Ar = ±2 (±14)kcV for A'47r (A'Gtt). 

Magnetic field. The magnetic field inside the track- 
ing volume has several components: the main solenoidal 
field, fields from permanent magnets and an additional 
magnetization of the latter due to the main field. To 



understand the effect of the field on the track recon- 
struction, the solenoid field strength is varied by ±0.02% 
and the r nag netization of the permanent magnets by 
±20% 0, if. For the mass difference, the larger de- 
viations arising from the change in rescaled solenoid 
field and magnetization are added in quadrature and the 
sum is assigned as a systematic uncertainty associated 
with the magnetic field; the same is done for the decay 
width. The resuffs are Aa™ = ±13 (±19) kcV/c^ and 
Ar = ±19 (±ll)keV for A:47r (ATBtt). 

Distance scale. A further source of uncertainty for 
the momentum determination comes from the distance 
scale. The positions of the signal wires in the drift cham- 
ber are known with a precision of 40 fj,ui, corresponding 
to a relative precision of 0.01%. As an estimate of the 
uncertainty of the momentum due to the distance scale, 
a systematic error half the size of the uncertainty ob- 
tained from the 0.02% variation of the solenoid field is 
assigned. For the mass difference this yields a shift of 
±4 (±6)keV/c2 for A:47r (ATGTr); the width is shifted by 
±8 (±4)keV for A'47r {K6Tr). 

Drift Chamber hits. In the standard D^^ selec- 
tion no lower limit is set for the number of drift cham- 
ber hits. Requiring at least 20 hits per track, thereby 
excluding the low momentum pions from D*~^ decays, 
modifies Am{D+^)o by -11 (-15)keV/c2 and T{D+^) by 
-7 (-7) keV for A'47r {K6Tr). 

Angular dependence. For the reconstructed AT^ 
and Z?" masses from the test data samples (see Sec. lIVp . 
a sine-like dependence on the azimuthal angle (f> is ob- 
served. This effect was also observed in a previous 
BaBAR analysis and might be related to the internal align- 
ment of the DCH [3^. For a detailed study, the same (p- 
dependence is introduced into the signal MC samples by 
modifying the kaon and pion track momenta accordingly. 
Due to symmetry, the effect disappears when all </> angles 
are taken into account, but as a conservative estimate the 
amplitude of the sine-like shift on the reconstructed D^^ 
mass in MC (13 (14) keV/c^ for A:47r {K6tv)) is taken as 
a systematic error for Am(_D^j^)o. 



IX. RESULTS 

For the combination of the measurements, a Best Lin- 
car Unbiased Estimate (BLUE, (s^j) technique is used, 
where correlations between the systematic uncertainties 
are taken into account. Adding the nominal and A'" 
masses, 2010.25 MeV/c^ and 497.614 MeV/c^ (with their 
respective errors of 0.140 MeV/c^ and 0.024 MeV/c^ ^), 
the final value for the Df^ mass is 

miD+^) = 2535.08 ± 0.15 McV/c^. 
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Using a slightly different method for the combination of 
the individual results a value of 

TOp+i) = 2535.08 ± 0.01 ± 0.15 McV/c^ 

is obtained, where the first error denotes the statistical 
and the second the systematic uncertainty. The latter 
is dominated by the uncertainty of the D*~^ mass. The 
mass difference between the D'^-^ and the D*^ is 

m{Dj^) - m{D*+) = 524.83 ± 0.04MeV/c^ 

using the BLUE technique, and for the alternative com- 
bination method 

m{D+^) - m{D*+) = 524.83 ± 0.01 ± 0.04MeV/c^ 

which has a significantly smaller systematic uncertainty 
than the m{Dfj^) result. 

For the total decay width of the D'^i, combining the 
results from the two measurements in the same way as 
for the mass yields 

T{Df^) = 0.92 ± 0.05 MeV, 

using the BLUE technique, and for the alternative com- 
bination method 

r{D+^) 0.92 ± 0.03 ± 0.04 MeV. 

The corrections of —51 (—45) kcV for the KAtt {K6tt) 
decay mode, based on the systematic resolution studies 
(Sec. IVIII Ap . are applied prior to the combination pro- 
cess. 



X. SUMMARY 

In this paper, precision measurements of the mass and 
decay width of the charmed-strange meson I?si(2536)+ 
via the decay D'^^ D*^K'^ are presented. Two decay 
modes are analyzed, with the that originates from the 
D*^ decaying either through K^n^ or if "tt+tt^tt". 

The results include the first significant measurement 
of the total decay width of the D^^ . This width is deter- 
mined to be 

r(D+) = 0.92 ±0.03 ±0.04 MeV, 



compared to the 90% confidence level upper limit of 
2.3 MeV given in Ref. [2l|. The mass of the £),i(2536)+ 
is measured to be 



m{D+^) = 2535.08 ±0.01 ±0.15 MeV/c^ 



and the D^^^ — mass difference to be 



m{D+^) - m{D*+) = 524.83 ± 0.01 ± 0.04MeV/c^ 



The result for the D+ - D*+ mass difference represents 
a significant improvement compared to the current world 
average of 525.04 ± 0.22 MeV/c^ [2l|. 

Based on a decay angle analysis, the = 1+ assign- 
ment for the D^^ meson is confirmed. 
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